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Quark-gluon vertex in arbitrary kinematics* 
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We compute the quark-gluon vertex in quenched lattice QCD, in the Landau gauge using an off-shell mean- 
field C(a)-improved fermion action. The complete vertex is computed in two specific kinematical limits, while the 
Dirac-vector part is computed for arbitrary kinematics. We find a nontrivial and rich tensor structure, including 
a substantial infrared enhancement of the interaction strength regardless of kinematics. 

1. INTRODUCTION 



Over the past few years, substantial progress 
has been made in our understanding of the non- 
perturbative correlation functions (propagators 
and vertices) of the fundamental fields of QCD 
and their relation to the phenomena of colour con- 
finement and dynamical chiral symmetry break- 
ing. It has recently become evident that at least 
in Landau gauge, a detailed knowledge of the 
structure of the quark-gluon vertex is essential 
for an understanding of the dynamics of quark 
confinement and chiral symmetry breaking, which 
is encoded in the quark Dyson-Schwinger equa- 
tion (DSE), relating the quark propagator S{p) to 
the gluon propagator and the quark-gluon vertex 
T/i (p, q) , where p and q are quark and gluon mo- 
menta respectively. 

The overall shape of the gluon propagator is 
now quite well known, both from lattice QCD 
|1I2| and from studies of the coupled ghost-gluon 
Dyson-Schwinger equations. It is now clear that 
if this is fed into the quark DSE together with 
a bare or QED-like vertex, the resulting quark 
propagator will not exhibit a sufficient degree of 
chiral symmetry breaking. Several of the con- 
tributions at this conference have addressed this 
issue, in various ways |8l4l5l()j . 

The quark-gluon vertex is related to the ghost 
sector through the Slavnov-Taylor identity (STI), 



g^r^(p,g) -G(g2) il-B{q,p + q))S-Hp) 



S-\p + q){l-B{q,p + q)) 



(1) 
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where G{q^) is the ghost renormalisation function 
and B°-{q, k) is the ghost-quark scattering kernel. 
In particular, if the ghost propagator is infrared 
enhanced, as both lattice 7 and DSE studies |HI 
indicate, the vertex will also be so. This provides 
for a consistent picture of confinement and chi- 
ral symmetry breaking at the level of the Green's 
functions of Landau-gauge QCD, where the same 
infrared enhancement that is responsible for con- 
finement of gluons, provides the necessary inter- 
action strength to give rise to dynamical chiral 
symmetry breaking in the quark sector. 

Confinement of quarks is still not understood 
in this picture, however. If the effective interac- 
tion between a quark and an antiquark by way 
of exchange of a nonperturbative gluon is to give 
rise to a linearly confining potential, the quark- 
gluon vertex must contain an infrared enhance- 
ment over and above that contained in the ghost 
self-energy. In the STI, this would be encoded 
in a non-trivial ghost-quark scattering kernel. 
Confirming or refuting this picture is a major 
challenge for current lattice and DSE studies of 
Landau-gauge QCD. 

Another area where the quark-gluon vertex 
may be of interest is that of effective charges. Al- 
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though 'the running couphng' is not a meaningful 
concept beyond perturbation theory, since there 
is no known way of nonperturbatively connecting 
two different 'schemes', process-dependent effec- 
tive charges may be defined non-perturbatively 
and be phenomenologically useful. The interac- 
tion between quarks and gluons may be a natural 
starting point for many of the physically interest- 
ing processes. 

Here we will present results of a lattice investi- 
gation into the quark-gluon vertex W . This con- 
sists of a determination of the full structure of the 
vertex at two particular kincmatical points (the 
soft gluon point where the gluon has zero momen- 
tum, and the quark reflection point where the in- 
coming and outgoing quark momenta are equal 
and opposite), as well as a determination of the 
dominant, vector part of the vertex in arbitrary 
kinematics. 

2. FORMALISM 

We denote the outgoing quark momentum p 
and the outgoing gluon momentum q. The in- 
coming quark momentum is k — p + q. In the 
continuum, the quark-gluon vertex can be de- 
composed into four components Li contributing 
to the Slavnov- Taylor identity and eight purely 
transverse components Tf. 



r^d?, q) 



i=l 



X,{p ,q ,k )L,^^,{p,q) 



(2) 



In euclidean space, the components Li and Ti are 
given by |3] 

=7m i2,M PP^ (3) 

Ti,M = - i^M 72,M =- 

Ts.M = - i'^M-g-' 76,M HlPht^- #-Pm (4) 
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where = p^+k^, = {pq)k^-{kq)p^. In Lan- 
dau gauge, for q ^ 0, we are only able to compute 
the transverse projection of the vertex, Tj^{p, q) = 
Pf^,yiq)^i'ip,q), where P^_,^iq) = S^^, - q^qvjq^ is 
the transverse projector. Since the vertex will 
always be coupled to a gluon propagator which 
contains the same projector, this is also the only 
combination that appears in any application. The 
four functions are projected onto the trans- 
verse li^^, giving rise to modified form factors 



A'l = Ai - g^Ta ; A2 = A2 — —r2 ; 
A3 = A3 —Tx ; A4 = A4 + T4 . 



(5) 



The lattice tensor structure is more complex, and 
(0) is only recovered in the continuum. The form 
factors also receive large contributions from lat- 
tice artefacts at tree level; the procedure we apply 
in correcting for these is described in 

In QED, the four form factors Aj are completely 
determined by the fermion propagator S'~^(p) = 
i i>A{p'^) -)- 

\r{p\q\k^) = \[A{f)-rA[e))- (6) 
A2(p,g,fc) = -- p ; (7) 



A3(p^g^fc2)^ 



,2 



p. 



M{p\q\k')^Q. 



(8) 
(9) 



By comparing our results with these forms we can 
get an idea of the importance of the nonabelian 
contributions to the STI, and in particular the 
quark-ghost scattering kernel. 

3. RESULTS 

We have analysed 495 configurations on a 16"^ x 
48 lattice at /3 = 6.0, using a mean-field improved 
SW action with a quark mass m « 115 MeV. 
This is part of the UKQCD data set described in 
|11| : further details can also be found in (9,. A 
second quark mass m w 60 MeV has also been 
studied, but as the mass dependence was found 
to be almost negligible , we do not show those 
results here. 
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3.1. Soft gluon point q = 

At g = the vertex reduces to 

r^(p,o) =Xi{p'h^ 



(10) 



where for brevity we write Ai(p^,0,p^) = K{p^)- 
In this specific kinematics, the QED expressions 
®-lIHll become^ 



(11) 
(12) 

(13) 



A?^°(p^)=A(p2); 

In Fig. n we show the dimensionless quantities 
Ai, 4p^A2 and 2pA3 as a function of momentum p. 
We also show the abelian forms (|11|I - H13|I which 
have been obtained from fitting lattice data for 
the quark propagator |1(JI12| . All these form fac- 
tors have been renormalised at 2 GeV, requiring 
Ai(4 GeV^0,4 GeV^) = 1. 

We find that while A3 is quite close to its 
abelian form, both Ai and A2 are significantly en- 
hanced. Since the ghost self-energy would con- 
tribute the same prefactor (in this kinematics, a 
constant) to all three form factors compared to 
the abelian form, this points to a nontrivial struc- 
ture in the quark-ghost scattering kernel. How- 
ever, the singular nature of the soft giuon point 
along with our small lattice volume make it diffi- 
cult to draw firmer conclusions. 

3.2. Quark reflection point p = ~k 

At the quark refiection point p = —k,q — —2p 
only A'j^ and T5 survive, and the projected vertex 
is 



r^(-<z/2,g) =X[{q'){^,- 



(14) 



where in this section we write 

{A'l,r5}(gV4,9^gV4) = {K.n}{q^)- The di- 
mensionless quantities A']^(q^), (7x5(9^) are shown 

■^In IIUI the expressions for A2 and A3 had the wrong sign. 
We are grateful to Craig Roberts for bringing these errors 
to our attention. In the same paper, the lattice data for 
A3 and for T5 at the quark refiection point also had the 
wrong sign. 




Figure 1. The renormalised, dimensionless form 
factors Ai, 4p^A2 and 2pA3 at the soft gluon point, 
as a function of p, for m — 115 MeV. Also shown 
are the corresponding abelian (Ball-Chiu) forms 
(|11|I - H13|I . derived from the quark propagator. 



in Fig. 12 These form factors have been renor- 
malised requiring A'l(lGeV^4GeV^lGeV^) = 1. 
A^ shows the same qualitative behavour as Ai at 
the soft gluon point, with a quite strong infrared 
enhancement. We see that rs, which has rarely 
if ever been included in vertex models used in 
DSE studies, is quite sizeable, indeed comparable 
in magnitude to the dominant component X'^. It 
will be interesting to study the effect of including 
this part of the vertex in future DSE studies. 

3.3. General kinematics 

The general lattice tensor structure, even for 
the Dirac-vector part of the vertex alone, is very 
complicated and makes a full determination of the 
vertex very difficult with this lattice action. How- 
ever, in the special case where we choose both the 
quark and gluon momentum vectors to be 'per- 
pendicular' to the vertex component, i.e. if we 
compute r^(p, q) withp^ = = 0, this structure 
simplifies considerably. There is no loss of gen- 
erality provided rotational symmetry is restored 
in the continuum. Here, we will only study the 
leading, vector part of the vertex, but the other 
components may also be determined in principle. 
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Figure 2. The renormalised, dimensionless form 
factors X'l and grs at the quark reflection point, as 
a function of the gluon momentum g, for m = 115 
MeV. 



Figure 3. The unrenormalised form factor X'l in 
the quark-symmetric kinematics = (up- 
per surface), as a function of quark momentum 
p (long axis) and gluon momentum q (short axis) 
in units of GeV. Statistical uncertainties are il- 
lustrated by the lower surface. 



In continuum notation, we compute 
1 



1 



5 

q 

{pq)kf, - (kq)pf, (p^ + kf^jX'^ 

-[k^-p^-ikl-pl)]r, 
X[ - [k^ ~ p')Te = X" . 



(15) 



(16) 



Of particular interest is the quark-symmetric 
limit, where the two quark momenta are equal in 
magnitude, p^ = k^. In this case, rg is also elim- 
inated, i.e. X'{{p'^,q^,p^) = A'i(p2,g2,p2). Note 
that both the soft gluon and the quark reflec- 
tion kinematics discussed previously are specific 
instances of this more general case. The details 
of the lattice implementation of this, including 
the tree-level correction, will be described else- 
where ■ In Fig- 0] we show A'^ as a func- 
tion of the two remaining independent momen- 
tum invariants. The data become quite noisy as 
q is increased, and also exhibit some 'spikes' and 
'troughs' which at present we assume to be nu- 
merical noise and lattice artefacts. 

By interpolating the points in Fig. 13 we may 
reach the totally symmetric point where p^ = 



k"^ = q^. This kinematics has a history of being 
used to define a momentum subtraction (MOM) 
scheme jjl]. We show our results in Fig.^J Again 
we find a strong infrared enhancement. 

Finally, figs. [3 and show A" in general kine- 
matics, for four different fixed values of q, as a 
function of the two quark momenta p and k. We 
expect all form factors to be symmetric in and 
k^ [tq on its own is antisymmetric, but is multi- 
plied by — fc^), and this is also what the figures 
show, within errors. The broadening of the data 
surface as q grows is simply a reflection of the 
increase in available phase space. 

The same qualitiative features as were found 
in the more specific kinematics, are reproduced 
here. At low q, we see a clear infrared enhance- 
ment, which disappears as q grows, reflecting the 
fact that at high momentum scales, only the log- 
arithmic behaviour (which is too weak to be seen 
in these data) remains. At the same time, the 
level of the surface sinks, which reflects the in- 
frared enhancement of A" also as a function of 
gluon momentum. 
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Figure 4. The unrenormalised form factor X[ at 
the totahy symmetric kinematics = = q^, 
as a function of the momentum p. 



4. DISCUSSION AND OUTLOOK 

We have determined the complete tensor struc- 
ture of the quark-gluon vertex at two kincmatical 
points, as well as the leading component in arbi- 
trary kinematics. At the soft gluon point, wc have 
observed significant and non-uniform deviations 
from the abelian form which has previously been 
the basis for DSE studies. At the quark reflection 
point, we find a significant contribution from the 
'chromomagnetic' form factor rs, which has pre- 
viously been ignored. In general kinematics, we 
find an infrared enhancement in all momentum 
directions; we hope to be able to quantify this en- 
hancement more clearly by fitting the data in the 
infrared region to functional forms in the three 
momentum variables (p^, /c^, g^). 

It is interesting to compare these results with 
recent calculations based on nonperturbative ex- 
tensions of the one-loop vertex |15l5j . Both these 
studies agree very well with our results for Ai and 
A3, while finding substantially lower values for 
A2. Since all these calculations must be consid- 
ered preliminary, too much emphasis should not 
be placed on this. It is worth noting that A2 is 
inherently noisy, as it mixes with Ai , which must 
be subtracted in order to obtain the data shown 
in Fig. n 




Figure 5. The unrenormalised form factor A'/ 
for gluon momentum q = 0.555 GeV (top) and 
q — 0.873 GeV (bottom), as a function of quark 
momenta p and k. The lower surfaces denote the 
statistical uncertainties. 



These results have been obtained on a rather 
small lattice, and with a discretisation that gives 
rise to quite large tree-level lattice artefacts which 
must be corrected for. We therefore expect sys- 
tematic errors to be large for large momenta. To 
obtain more reliable results, and to extend this 
study to the full vertex structure at all kine- 
matics, it would be desirable to employ an ac- 
tion which is known to have smaller and more 
tractable tree-level artefacts. The Asqtad action 
has been employed successfully in computing the 
quark propagator j and unlike the SW action. 
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Figure 6. As Fig. |S1 but for gluon momentum 
q = 1.193 GeV and 2.321 GeV. 



only Ai and possibly A2 are non-zero at tree level, 
so tree-level correction will not be needed for the 
remaining form factors. This action is also com- 
putationally cheap, making large lattice volumes 
feasible. Another possibility is to use overlap 
fermions, which have the advantage of retaining 
an exact chiral symmetry, which protects all the 
odd Dirac components of the vertex at tree level. 
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